New description of four-body breakup reaction 
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We present a novel method of smoothing discrete breakup cross sections calculated by the method of 
continuum-discretized coupled channels. The smoothing method based on the complex scaling method is tested 
with success for a ^**Ni(d, pn) reaction at 80 MeV as an example of three-body breakup reactions and applied to 
a ^^C(''He, nn"*He) reaction at 229.8 MeV as an example of four-body breakup reactions. Fast convergence of 
the breakup cross section with respect to extending the model space is confirmed. The method is also applied to 
^^C(''He, Tin^He) and ^°*Pb(®He, nn^He) reactions at 240 MeVM and compared with the experimental data. 
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Exploring unstable nuclei far from the stable line is one of 
the most important subjects in nuclear physics. The unstable 
nuclei have exotic properties such as the halo structure 
[st] and the island of inversion As a feature of reactions 
induced by unstable nuclei, the projectile easily breaks up 
into its constituents. One of the most reliable methods for 
treating the projectile breakup processes over a wide range of 
incident energies is the continuum-discretized coupled chan- 
nels (CDCC) method iHISl. In CDCC, the scattering wave 
function of the total system is expanded with a finite number 
of bound and discretized continuum states of the projectile. 
The space spanned by these states is called the model space. 
The S'-matrix elements calculated with CDCC converge as the 
model space is extended 10,11!]. The converged CDCC solution 
is the unperturbed solution of the distorted Faddeev equations, 
and corrections to the solution are negligible within the spatial 
region in which the breakup processes take place Igl fioll . 

For scattering of a two-body projectile, the continuum 
states are classified by linear and angular momenta, k and /, 
respectively, between the two constituents. In CDCC, these 
momenta are taken up to upper limits, the k continuum is di- 
vided into small bins and the continuum states in each bin 
are averaged into a single state. This discretization proce- 
dure is called the average (Av) method. The Av method has 
been widely used, but its application has been hmited to three- 
body breakup reactions as we will show. An alternative to 
the Av method is the pseudostate (PS) method lfTll - [l6ll . in 
which the continuum states {ip'^^^k)} are replaced by pseu- 
dostates {$„} obtained by diagonalizing the internal Hamil- 
tonian of the projectile in a space spanned by i^-type basis 
functions. One can adopt the transformed harmonic oscilla- 
tor (THO) m or the Gaussian |[ll[ll as the i^.type basis 
functions. The validity of the PS method was confirmed for 
scattering of two-body projectiles by the agreement between 
CDCC solutions calculated with the two discretization meth- 
ods lE 
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For scattering of ®He as a typical example of four-body 
breakup reactions, CDCC with the PS method based on Gaus- 
sian lfT4l [Tsll or THO lfl6ll basis functions was successful in 
describing the elastic scattering at not only high energies but 
also low energies near the Coulomb barrier Thus, the back- 
coupling effect of four-body breakup processes on the elastic 
scattering is well described by the PS method. 

For the three-body projectile, continuum wave functions 
-0'"' (p, fc) are classified by momenta (p, k) conjugate to two 
internal coordinates (p,r) of the three-body system. The 
breakup 5-matrix elements calculated with CDCC, 5„, are 
discrete in p and fc, although the exact ones S{p, k) are 
continuous. Thus, one needs a way of smoothing Sn- In 
principle this is possible by calculating the smoothing factor 
(■(/''^■'(p, fc)|"I>n) l,12il , but in practice it is not easy, because 
evaluating ^^^^{p, k) for many combinations of p, k and n 
is quite time consuming. Recently, we proposed two methods 
of obtaining the smoothing factor; one with a direct numeri- 
cal integration lUTIl and another [TSfl with the complex scaling 
method (CSM) 111911 . However, these require complicated nu- 
merical calculations, so that the convergence of the smoothed 
breakup cross section with respect to increasing the model 
space is not sufficient particularly for the differential breakup 
cross section as a function of the excitation energy e of projec- 
tile, da/de. If the Av method is applied with small momentum 
or energy bins, one can obtain da/de without evaluating the 
smoothing factor. This was done with the hyperradial contin- 
uum wave function ll20ll . but the convergence of the CDCC 
solutions has not been obtained yet for da/ds 1I21I1 . A way 
of circumventing these difficulties is to construct a method of 
obtaining the smoothing factor without evaluating t/''^' (p, k). 

In this Rapid Communication, using CSM and CDCC, we 
propose a practical method of obtaining da/de as a contin- 
uous function of e without evaluating -(/^'^'(p, fe). CSM is a 
powerful tool for obtaining many-body resonance and weakly 
bound states |f2^. Recently, it was applied to the electromag- 
netic transition of the core H- nucleon H- nucleon system, such 
as ^He and ^^Li from the ground state to the continuum state 
with e ll23l l24ll . CSM is applicable not only for resonances 
but also for continuum states, so that the transition strength 
is obtained as a continuous function of e. The smoothing 
method proposed here is an alternative to the direct calcu- 
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lation of the smoothing factor fc)| j>„) with the hy- 

perradial continuum wave function 11251 [2611 . The former is 
considered to be more practical than the latter, because the 
former does not require one to evaluate the continuum states 
k) for many combinations of p and k. The validity 
of the new method is tested for a three-body breakup reac- 
tion, ^*Ni((i, pn) at 80 MeV, in which the "exact" breakup 
cross section is obtainable by calculating the smoothing fac- 
tor with the direct numerical integration. The new method 
is applied to the ^^CC'^He, nn'^He) reaction at 229.8 MeV. A 
merit of the present smoothing method is that one can see fast 
convergence of the calculated breakup cross section with re- 
spect to extending the model space. The method is also ap- 
pHed to i2c(6He, nn^He) and ^°^Ph(^He, nn^He) reactions 
at 240 MeVM and compared with the experimental data. In 
principle, this method is applicable not only for four-body 
breakup reactions but also for many-body breakup reactions. 

We consider scattering of a projectile B from a target A, 
in which B is composed of three constituents (x = b, c and 
d). The scattering is described by the four-body Schrodinger 
equation 

(i7-£;tot)|«'^+^) =0 (1) 

with the outgoing boundary condition, where the total energy 
E'tot satisfies i^tot = E^^+Eq for the corresponding incident 
energy in the center of mass of the B + A system and 
the ground-state energy Eq of B. The total Hamiltonian H of 
this system is defined by 



H ^Kr + U + Hb (2) 

with 

U = Ub + Uc + Ud + 14^°"' + + , (3) 

Hb = Ky + Kr + Vbc + Vcd + Vdb , (4) 



where Hb is an internal Hamiltonian of B. The relative co- 
ordinate between B and A is denoted by R, and the internal 
coordinates of B are denoted by a set of Jacobi coordinates 
^ = (y, r). Momenta conjugate to coordinates R and (y, r) 
are represented by P and {p, k), respectively. The kinetic 
energy operator associated with R (^) is represented by Kr 
(K^), Vxx' is a nuclear plus Coulomb interaction between x 
and x', and Ux and V^°^^ are nuclear and Coulomb potentials 
between x and A, respectively. 

In CDCC with the pseudostate discretization method, the 
scattering is assumed to take place in a model space lfl2l - [l5ll : 

7' = ^|$„)(<I>„|, (5) 

n 

where (&„ is an ??th eigenstate obtained by diagonalizing Hb 
with L^-type basis functions. For simplicity, B is assumed 
to have only one bound state $0. The four-body Schrodinger 
equation is then solved in the model space: 

ViH - E,,,)P\¥+^^^) = 0. (6) 



The model space assumption has already been justified by the 
fact that calculated elastic and breakup cross sections con- 
verge with respect to extending the model space EMI. 

The exact T-matrix element to a breakup state with (p, k) 
can be described by 

T,{p,k,P) = {4-\p,k)xi-\P)\U-V§°^'\¥+^), (7) 

where Vg °"' is a sum of Coulomb interactions between B and 
A but the arguments are replaced by R: 

V§°''\R) = 14'^°"' (i?) + + (8) 

The final-state wave functions, jV's \p,k)) and \xi \P)), 
with the incoming boundary condition are defined by 

[Tn + V§°^\R) - (Stot - e)] Ixi"' (P)) - 0, (9) 
{HB~e)\4-\p,k)) = 0, (10) 

where Sft - e ^ {hP)^/{2nR) and e = {hp)^/{2ny) + 
(hk)'^/ (2/ir) for reduced masses fiR and p^ of coordinates R 
and ^, respectively. Inserting the approximately complete set 
Eq. Q into Eq. (jT]), we can find lfl2l - [l5ll that the T-matrix 
element is well approximated by 

r,(p,fe,P)«^(V'^)(p,fc)|<i>„)r„ (11) 

with the CDCC T-matrix element 

to an jith discrete breakup state $„ with an eigenenergy e„. 

Here Eq. ( fTTT ) is derived by replacing P by P„ in Xe \P)- 
Tn is obtainable by CDCC, but it is quite hard to calculate the 
smoothing factor {il>i ■* (p, fc) |$„) directly with either numer- 
ical integration iflTIl or CSM ifisll . Hence, we propose a new 
way of obtaining the differential cross section with respect to 
e without calculating the smoothing factor 

Using Eq. ( fTTT i. one can rewrite the differential cross sec- 
tion as 

= ldp'dk'd{s-s')\TAp',k',P')f 

« -nie,np) (13) 

TT 

with the generalized response function 

7^(e,^^p) =Im [ J2 T*{'fn\G^-'^\^n')TnA, (14) 

where G^~^ = lim (e - Hb - ir])^^- In Eq- ( fT4] l. there 

is no smoothing factor, as expected. Furthermore, the propa- 
gator G^^^ operates only on spatially damping functions $„, 
so that the calculation of ($„|G''^-' |<E>,i') becomes feasible, as 
we will show. 
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CSM is now applied to evaluating ($„|G(~) |$„'). The 
scaling transformation operator C{d) and its inverse are de- 
fined by 

{y,r\Cie)\f) = e^^Vlye^^r-e'"), (15) 

{f\C-\9)\y,r) - [e-^^'fiye-''^,re-^')Y . (16) 

Using the operators, one can get 

($„|G(-'|$„0 = ($„|C^i(0)G^-^C(0)|$„.), (17) 



where 



G 



(-) 



lim — 

v-^+o e ■ 



m 



11] 



(18) 



with Hfj = C{e)HBC~^{0). When < < 0, the scaled 
propagator (^|Gg~'' |^') is a damping function of ^ and It 
should be noted that although the scaling angle in general cal- 
culations with CSM has been taken as positive, the angle in 
the present situation becomes negative since G^^' has the in- 
coming boundary condition. Hence, it can be expanded with 
L^-type basis functions with high accuracy: 



G. 



(-) 



(19) 



where (j)^ is an ith eigenstate obtained by diagonalizing Hg in 
a model space spanned by L^-type basis functions: 



(20) 



By virtue of CSM, thus, we do not need to evaluate 
the exact three-body continuum state i/'e \k,p) to obtain 

(<1>„|G(-)|$„0- 

Inserting Eq. iT% into Eq. ( fT4b by means of Eq. (T% leads 
to a useful equation: 



dedili 



1 V- 

-Im > 



with 



n' 



(21) 

(22) 
(23) 



The principal result of this Rapid Communication is that 
C{9) and C~^{9) operate only on the spatially damping func- 
tion $„. This makes the calculation of (0^|G(6')| <!>„') and 
{^n\G~^ {0)\(j)f) feasible and makes possible the convergence 
of d^a/dedflp with respect to extending the model space, as 
we will show. In other words, C{9) and C^^{d) are not al- 
lowed to act on a nondamping function such as the plane wave, 
since the scaled function diverges asymptotically in this case. 

We now test the validity of Eq. d^TT i for a three-body 
breakup reaction in which the "exact" breakup T-matrix ele- 
mentT(fc,P) = X)„(V'*~-'('s)l'J'n)^n is obtainable by taking 



the overlap (■!/''"■' (fe) |$„) directly with numerical integration. 
This approach is referred to here as the "smoothing" factor 
method ||T3l. As an example, we consider a ^*Ni((i, pn) re- 
action at 80 MeV; see Ref. Ill] for the details of the CDCC 
calculation. Figure[T]shows the differential breakup cross sec- 
tion da/de in which the double differential cross section of 
Eq. ( l2n i is integrated over the solid angle fip of momen- 
tum P. The new method (open circles) yields the same result 
as the smoothing factor method (solid line). Thus, the new 
method is confirmed to be valid. 




10 
£ [MeV] 



FIG. 1; Calculated breakup cross section for d + ^^Ni scattering at 80 
MeV. The solid line shows the result of the smoothing factor method; 
open circles show the result using the new method. 



Next, the new method is applied to ^^C(^He, nn^Yie.) scat- 
tering at 229.8 MeV as an example of a four-body breakup re- 
action; see Ref. fl6ll for optical potentials between n and ^^C 
and between a and ^^C. As interactions Vnn and Vna in Hb, 
we take the so-called GPT iH and KKNN 128:] potentials, 
respectively. These potentials with a Gaussian form repro- 
duce well data of low-energy nucleon-nucleon and nucleon- 
^He scattering, respectively. The particle exchange between 
valence neutrons and neutrons in "'He is treated approximately 
with the orthogonality condition model ll29ll . 

For the diagonalization of Hb and Hg, we adopt the Gaus- 
sian expansion method (GEM) f3(|. In GEM, the state of the 
^He + n + n system is described by a superposition of three 
channels, each channel with a different set of Jacobi coordi- 
nates, (y^,rc). For each channel (c), the radial parts of the 
internal wave functions involving and Tc are expanded by 
a finite number of Gaussian basis functions 

-(yc/i/j) 



(24) 



respectively. Here \{£) is the angular momentum of ydvc), 
and the range parameters are taken to lie in geometric pro- 
gression: 



(25) 
(26) 



The parameters depend on c, but we omitted the dependence 
in Eqs. ( |25] ) and ( |26] | for simplicity; see Ref. |[l4ll for the 
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details of the diagonalization and the definition of Jacobi co- 
ordinates. 



TABLE I: Gaussian range parameters. 



Set 


c 


J max 


Vi (fm) y 


max (fm) 


^max 


ri (fm) 


rmax (fm) 


I 


3 

1, 2 


10 
10 


0.1 
0.5 


10.0 
10.0 


10 
10 


0.5 
0.5 


10.0 
10.0 


II 


3 

1, 2 


15 
15 


0.1 
0.5 


20.0 
20.0 


15 
15 


0.5 
0.5 


20.0 
20.0 


m 


3 

1, 2 


20 
20 


0.1 
0.5 


50.0 
50.0 


20 
20 


0.5 
0.5 


50.0 
50.0 




£ [MeV] 



In order to confirm the convergence of the breakup cross 
section with respect to extending the model space, we pre- 
pared three sets of basis functions shown in Table |T] For 0"*" 
and 1~ states, maximum internal angular momenta are set to 
i = \ = 1. For 2+ states, they are = A = 1 for c = 1 and 2, 
and £ = A = 2 for c = 3. 



FIG. 3: Convergence of the breakup cross sections to 2 continuum. 
Definition of lines is the same as in Fig.|2] 
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FIG. 2: Convergence of the breakup cross sections to (a) 0^ con- 
tinuum and (b) 1~ continuum. Dashed lines, solid lines, and open 
circles correspond to results of sets I, II, and III, respectively. The 
dotted line in (b) shows the results when Coulomb breakup processes 
are switched off. 



Figure[2]shows the breakup cross sections da/de to 0+ and 
1~ continua separately, and Fig. [3] shows the cross section 
to the 2+ continuum. For all the cross sections, sets II and 
III yield the same result, but the result of set I is somewhat 
different from those of sets II and III. Thus, the convergence 
with respect to increasing the model space is obtained with set 
II. Figure |4] shows the dependence of the net breakup cross 
section to 0+, 1^, and 2+ continua around the 2+ resonance 
peak. The net breakup cross section converges at 6* = —14° 
when 6 decreases from —6° to —18°. 

The present calculation includes Coulomb breakup pro- 
cesses that were neglected in the previous paper ITsll . The 
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FIG. 4: Dependence of the net breakup cross section to 0"'", 1^, and 
2^ continua on 9. 



effect is significant for the 1^ continuum, as shown by the 
dotted line in Fig. I2b). It enhances the breakup cross sec- 
tion by a factor of 2 from the result without Coulomb breakup 
processes. The effect is negligible for 0+ and 2+ continua, 
although the latter is a dominant component of the breakup 
cross section. Consequently, the Coulomb breakup effect 
in the present reaction system is not significant for either 
the breakup reaction or elastic scattering. For heavy targets. 
Coulomb breakup processes dominate breakup reactions and 
Coulomb breakup reactions are a useful tool for investigating 
properties of halo nuclei. The new method proposed here can 
treat both nuclear and Coulomb breakup processes and then 
be used to analyze their interference in the same framework. 

In Fig. |5] the breakup cross section da/de calculated 
by the present method is compared with the experimen- 
tal data for ^He + and *^He + ^o^pb reactions at 240 
MeVM ifsill . These data have aheady been analyzed by four- 
body distorted-wave Born approximation (DWBA) Il32ll and 
the eikonal approximation [33]. In the present calculation, we 
take the same potentials as in Ref. ll33ll for n-^"^Pb and a- 
^°^Pb systems. The optical potential for a n-^^C system is 
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e [MeV] 

FIG. 5: Comparison of the breakup cross section calculated by 
CDCC (solid line) with experimental data for (a) ®He + ^^C scat- 
tering at 240 MeV/A and (b) "He + ^'*Pb scattering at 240 MeVM. 
Dot-dashed, dotted, and dashed lines show contributions of 0^, 1^, 
and 2+ breakup, respectively. The experimental data are taken from 
Ref. (H. 



taken from the global nucleon-nucleus potential ll34ll . while 
the optical potential for a-^^C system is constructed from the 
+ potential at 200 MeV/A 13511 by changing the radius 
parameter from ^^C to a. Nuclear breakup is dominant for 
^He + ^^C scattering at 240 MeV/A, while Coulomb breakup 
to the 1^ continuum is dominant for ^He + ^osp^ scattering. 



For a ^^C target, the present theoretical result is consistent 
with the experimental data except for the peak of the 2+ res- 
onance around e = 1 MeV. This overestimation is also seen 
in the results of four-body DWBA, and the problem is partly 
solved by considering the experimental energy resolution. For 
208p|^ target, the present method underestimates the experi- 
mental data at £ > 2 MeV. A possible origin of this underesti- 
mation is that the inelastic breakup reactions are not included 
in the present calculation. As mentioned in Ref. ll32ll . the in- 
elastic breakup effect is not negligible, and the elastic breakup 
cross section calculated with four-body DWBA also underes- 
timates the data. 

In summary, we have proposed a practical method of cal- 
culating the differential breakup cross section as a continu- 
ous function of the excitation energy of a projectile, by com- 
bining CDCC and CSM. This method does not require one 
to calculate the continuum wave functions of the projectile. 
All we have to do is just diagonalize the projectile Hamilto- 
nian and the scaled Hamiltonian with L^-type basis functions. 
In the present formalism, the scaling operator C{9) operates 
only on spatially damping functions and hence the differential 
breakup cross section converges quickly as the model space is 
extended. The method is successful in reproducing the data 
on ^He + and ^He + ^ospb reactions at 240 MeVM. 
In principle, the present formalism is applicable for many- 
body breakup reaction, if the diagonalization of the projectile 
Hamiltonian and the scaled Hamiltonian is feasible. 
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